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Abstract 

The bound state of constituent quarks forming a Qqq composite baryon 
is investigated in a QCD-inspired effective light-front model. The light-front 
Faddeev equations are derived and solved numerically. The masses of the spin 
1/2 low-lying states of the nucleon, A , A+ and Ajj are found and compared 
to the experimental data. The data is qualitatively described with a flavor 
independent effective interaction. 
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I. INTRODUCTION 



One major task in strong interaction physics is the calculation of the wave function and 
the spectrum of the hadrons from Quantum Chromodynamics (QCD) HJ. Phenomenological 
dynamical models that retain the low-energy physics of QCD, and relate different observables 
are still of interest as long as they allow to expose the nonperturbative properties of QCD. 
One possibility to model the quark dynamics within a relativistic framework is to use the 
light-front dynamics in a truncated Fock-space, which yields a wave-function covariant under 
kinematical boosts In general, the light-front Fock-state truncation is stable under 

kinematical boost transformations @|. 

A light-front QCD-inspired model was recently applied to the pion and other mesons 
. It was able to describe the pion structure as well as the masses of the vector and 
pseudo-scalar mesons. This effective model, the ||-model ||, has two components in the 
interaction, a contact term and a Coulomb-type potential. The contact term is essential to 
collapse the constituent quark-antiquark system to form the pion, while the vector meson 
is dominated by the Coulomb-type potential. In this model the vector meson corresponds 
to a weakly bound system of constituent quarks. The model has no confinement, and the 
spin does not play a dynamical role besides justifying the contact term from the hyperfine 
interaction. However, the contact term was able to embrace the physical scale brought by 
the pion mass and from that the masses of the other pseudo-scalar and vector mesons were 
calculated ||. The flavor independence of the interaction was assumed, and the model 
reproduced fairly well the data, despite its simplicity. 

Here we report a nonperturbative calculation of the flavor dependence of the masses of 
some baryons, extending the concepts coming from the ||-model applied previously only to 
mesons |5|,^] . In this work, we show that the flavor independence of the effective interaction 
still holds for baryons, by studying the spin 1/2 low- lying states of the nucleon, A , A+ and 
A°. To get insight on the complex three-quark relativistic dynamics in this first study of 
the Qqq system within the light-front framework, we use only the contact interaction, which 
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brings the physical scale of the ground state of the nucleon, while the spin is averaged out. 
The effect of the long-range Coulomb-type interaction is effectively carried out by tunning 
the contact interaction to the nucleon mass. The mass of one of the constituent quarks (Q) 
will be varied while the bare strength of the effective contact interaction is kept constant. For 
each constituent mass the binding energy of the three-quark system is evaluated. Naturally, 
this calculation yields the binding energy of the constituent quarks as a function of the 
baryon ground state mass, because the binding energy and the ground state mass, depend 
only on the mass of the quark Q, with the other inputs kept unchanged. 

The binding energy for constituent quarks is a difficult concept to use together with quark 
confinement which is believed to exist in nature. It is fair to ask, if one has a meaningful 
model without confinement, how one could extract from the experimental data a quantity 
to be compared with the model binding energy. This key point will be addressed in detail 
in the text. 

Let us mention that, the three-body model with constituent quarks interacting in the 
light-front with a contact force |7| has been applied to the proton and described its mass, 
charge radius and electric form factor up to 2(GeV/c) 2 ||, although the spin was averaged 
out and before the advent of the j j-model Recently the same three-quark model was 

applied to study the dissolution of the nucleon at finite temperature and baryonic density 

The paper is organized as follows. In Sec. II, we derive the coupled integral equations 
for the Faddeev components of the vertex of the three-body light-front bound-state wave 
function 0, which generalizes the Weinberg- type integral equation found for a two-body 



bound state JTTJ. In Sec. Ill, we present the numerical results for the masses of the nucleon, 
A , A+ and A°, using the flavor independent contact interaction. We show as well how to 
relate to experimental data the results for the binding energies. Finally, in Sec. IV, we give 
our conclusions. 
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II. THREE-QUARK RELATIVISTIC MODEL 



The light-front is defined by x + = x° + x 3 = and the coordinates in this space- 
time hypersurface are given by x~ = x° + x 3 and ) MM- The coordinate 
x + is recognized as the time and k~ = k° — k 3 , the momentum canonically conjugated, 
corresponds to the light-front energy. The momentum coordinates k + and k±, are the 
kinematical momenta canonically conjugated to x~ and x±, respectively. 

A relativistic model for three-particles on the light-front for a pairwise contact inter- 
action, was derived from the three-body ladder Bethe-Salpeter equation for the Faddeev 
component of the vertex function, by eliminating the relative x + -time between the particles 
0. The projection of the covariant dynamics to the light-front hypersurface, is performed 
through the integration over the k~ momentum of the individual particles in the Bethe- 



Salpeter equation leading to a Weinberg-type equation |TTJ for three-particles. In the work 
of Ref. 0, in fact the kernel of the Faddeev equation in the light-front was derived in low- 
est order, and in principle corrections of higher order can be systematically constructed 



following recent discussions |l2-r4|. 

Let us sketch the derivation of the light-front Faddeev equations for a heavy-light-light 
three quark system (Qqq) from the four- dimensional Bethe-Salpeter equations. Two differ- 
ent spectator functions, which corresponds to the Faddeev components of the vertex, are 
possible. For the interacting pair being qq the spectator function is vq(q ,m ), function of the 
four-vector momentum of the quark Q. For the interacting pair Qq the spectator function 
is Vq^). The coupled Faddeev-Bethe-Salpeter equations in the ladder approximation are 
given by: 

d 4 k v q (gO 



Trig + 16)((Pb — q — k) 2 — m 2 q + it 
which is represented diagrammatically in figure 1, and 



v q (q") = ^r Qq (M 2 Qq ) J 



d 4 k 



v q (^) 



(k 2 — m 2 + ze)((Pe — q — k) 2 — itlq + is) 

v q m 

ik 2 — m 2 Q + i€)((Pb — q — k) 2 — m 2 + te) 
4 



(2) 



which is represented in figure 2. The baryon four- momentum is given by P#, the light 
and heavy quark masses are m q and uiq, respectively. The masses of the virtual two- 
quark subsystems are M qq = (Pg — q) 2 and Mq = (Pb — q) 2 due to the conservation of 
the total four-momentum. The two-quark scattering amplitudes r qq [M qq j and TQ q [Mq ) 




are the solutions of the Bethe-Salpeter equation in the ladder approximation for a contact 
interaction between the quark pairs, which are derived in detail in the appendix. 

The analytical integration over k~ is performed in Eqs. ([!]) and (§), using only the pole 
of the single quark propagator in the lowest half of the complex k~ plane. The pole is given 
by the on-energy- shell condition k~ n = (£r[ + m a)/k + , with a = q or Q when v q or 
vq (q 11 ) is integrated, respectively. The condition for a nonvanishing result of the integration 
in A;" is < k + < Pg — q + . The spectator functions appearing inside the integrations in 
Eqs. (0) and (0) depend only on the kinematical momentum, (k + ,k±), as long as k on is a 
function of the kinematical momentum. Thus, to close the light-front equations the external 
momentum is chosen on the fc~-shell. 

In the rest-frame of the baryon of mass M B , we write that v a (q + ,q±,q~ n ) = v a (q±,y), 
where for convenience the Bjorken momentum fraction y = q + /Mb is used. The Faddeev 
equations in the light-front written in terms of the kinematical momenta are given by: 






y 



X 



l-x-y 



(3) 




(4) 
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where x = k + /Mb- In the first equation of the coupled set, Eq. (0), Q is the spectator quark 
and qq is the interacting pair. The maximum value for k± is chosen to keep the mass squared 
of the qq or Qq subsystem real, i.e., M 2 > and Mq q > 0, respectively. These constraints 
in the spectator quark phase-space come through the theta functions in the integrations 



of Eqs. (§) and (|). For M 2 q > one has k± < k^{m q ) = J(l-x)(M B x - m 2 q ), and 



x > {m q /M B ) 2 . For M 2 q > one has k± < £^ ax (m Q ) = J(l - x)(M B x - m 2 Q ), and 
x > (wq/Mb) 2 . For equal particles, Eq.(||), reduces to the one derived in Ref. 0. 

Finally, the light-front baryon bound state wave-function of the Qqq system in the rest- 
frame is constructed from the Faddeev components of the vertex as in Ref. |§: 

v q fa, h±) + v q (x 2 , k 2 ±) + v Q (x 3 , k 3± 
^/xTx^(M 2 - M$) 



k 1± ; x 2 , k 2 ±) = — '— = ^ , 2 2 , , (5) 



where the free three-quark mass Mq is given by: 

M 2 = _1± q + _2± q_ + _6± Q _ /gN 

Xi x 2 x 3 

Each constituent quark has momentum fraction Xj and transverse momentum kj±(j = 1,3), 
satisfying X\ + x 2 + x 3 = 1 and ki±_ + k 2 ± + k 31 _ = 0. 



III. RESULTS 

The coupled integral equations (0) and (|J) for a relativistic system of three constituent 
quarks with a pairwise zero range interaction, are solved numerically. The physical inputs 
of the model are the constituent quark masses and the diquark bound state mass. The 
mass of the ground state baryon (M B ) and the binding energy, B B = 2m q + rriq — M B , are 
calculated. First in this section, before presenting the model results, we provide a qualitative 
discussion in order to attribute to the low- lying spin 1/2 baryons a binding energy from the 
experimental data. Then, we compare these data with the model calculations. 
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A. Qualitative analysis 



According to the effective QCD-inspired model calculations of Ref. IJ, the low-lying 
vector mesons are weakly bound systems of constituent quarks while the pseudo-scalars 
are more strongly bound. This justifies our supposition that the masses of the constituent 
quarks can be derived directly from the vector meson masses: 



m u 


= -M p = 
2 p 


0.384Ge1/ , 


m s 


= M K * - 


^M p = 0.508GeV , 


m c 


= Md* - 


-M p = 1.623GeV , 
2 p 


m b 


= Mb* - 


^M p = 4.941GeV , 



(7) 

where the values from Table I are used. Also one can check whether the values of the 
current quark masses obtained from the constituent ones, attain values compatible with the 



actual knowledge [IB]. It is reasonable to think that the constituent quark mass formation 
mechanism does not distinguish in detail the quark flavor and thus, the current quark mass 
of s, c and b are just given by the difference m c q rr = niq — m u (Q = s, c, b). The extracted 



values of the current masses are quite consistent with the experimental ones from Ref. |IB 



as one can verify in Table I. At least from this point of view is not unaceptable to define the 
constituent quark masses from the low-lying vector meson states. 

A remark should be added on how ambiguous the qualitative estimates are for the con- 
stituent and current quark masses. Errors arise when different meson masses are used as 
input, and from the current quark masses of the up and down quarks, which we have disre- 
garded and leads to an error of about 10 MeV. Moreover, another 10 MeV can be attributed 
to the degeneracy of the p and u mesons in our model. Therefore, we roughly estimate an 
error of 20 MeV in Table I for our values of current quark masses, and for the constituent 
quark masses from Eq. (]?]), as well. 

Below, we attribute values to the baryon binding energies, using the constituent quark 
masses from Eqs. (0) and the experimental values of the baryon masses jl6 ]: 
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r>exp 




- 3 -M - 
2 p 


- M 










= 2m u + m s 


71 /f 

— M^o = 


i\ /f 


2 p 


/I /f 

- M A o , 




n exp 


= 2m u + m c 


- M kt = 


M D * 


+ iM„ 






Tjexp 




= 2m u + mj 


-M A o = 


M B * 




-M A o . 


(8) 



The results are presented in Table II. 

In figure 3, the binding energies of the low-lying pseudo-scalar mesons [Bm\ defined as 
the difference between the vector and pseudo-scalar masses (see Table I) are plotted against 
the mass of the corresponding pseudo-scalar meson. Also, the values of the binding energies 
of the spin 1/2 baryons (N, A , A+ and A°) from Eq.(|8|) are shown as a function of the 
corresponding baryon mass. We observe a smooth trend in the plot of figure 3 for Bm as 
well as for as a function of the hadron ground state mass. The increase of the heavy 
quark mass produces the same qualitative behaviour irrespective of the nature of the hadron, 
being a meson or a baryon. 

The data for mesons shown in figure 3, was described by the effective QCD-inspired 
model once the hypothesis of the flavor independence of the interaction was adopted M. 
This is consistent with the fundamental QCD-theory in which the gluon does not recognize 
flavor but color |]J. Figure 3 suggests that for baryons, it is reasonable to assume, as a first 
guess, that the constituent quark interaction would be flavor independent. 

B. Model calculations 

The physical input of the light-front model defined by the coupled integral equations @ 
and (|4]) are the constituent quark masses and the diquark bound state mass. Up to this 
point the masses of the constituent quarks have been defined by Eq. (0). The value of the 
diquark mass has to be found. The diquark mass was fitted to the value of the proton mass, 
with a given constituent quark mass by the solution of Eq.(Q) with q = Q = u. We use 
m u = 0.386 GeV which together with the nucleon mass of 0.938 GeV implies in Md = 0.695 
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GeV. The slightly different m u in respect to the one found in Eq.s (0) is just to adequate 
to the value obtained in Ref. ||, where the above quark mass was used with a reasonable 
description of the proton charge radius and electric form factor below 2(GeV/c) 2 . 

We solve the coupled equations (0) and (f|) for fixed m u = 0.386 GeV and Ma =0.695 
GeV and different values of tuq. Thus, the binding energy for the spin 1/2 baryons A , A+ 
and A° are obtained by changing the value of trq (Q = s, c, b). Each value of tuq produces 
a ground state mass and binding energy. In figure 4, instead of showing the binding energy 
as a function of tuq, we plot it as a function of the ground state mass in a continuous 
curve and compare with the attributed experimental values. For the baryon mass above 
2.3 GeV, the bound Qqq system goes to the diquark threshold. This gives the saturation 
value of Q.QUGeV seen in figure 4. The model calculation with a flavor independent ef- 
fective interaction is able to reproduce the trend of the attributed experimental binding 
energies as a function of the mass of the baryon ground state. In view of the simplicity 
of the model, the agreement between the binding energies obtained theoretically with the 
attributed experimental values is quite reasonable. 

Our present results, although, in a simplified model generalizes the flavor independent 
effective interaction of the QCD-inspired |l-model to the context of the dynamics of con- 
stituent quarks forming the baryon. The very existence of the smooth pattern shown in 
figures 3 and 4, for the correlation between binding and masses, tells us that the dominant 
physics is the mass variation of the constituent quark and the spin effects should average 
out. Taking into account that spin degree of freedom is averaged out in the model and at 
the same time the reasonable description of the data found in figure 4, give us the confidence 
that the main physics related to quark mass variation is reasonable described by the flavor 
independent contact interaction. 
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IV. CONCLUSIONS 



We have studied the binding of the constituent quarks forming the the low-lying spin 
1/2 baryonic states of the nucleon, A , A+ and A°. We use a relativistic three-quark model 
of the baryon defined on the light-front, where the inputs were the constituent quark masses 
and the diquark mass in the light sector. The effective interaction was chosen of a contact 
form and spin was averaged out. The contact interaction includes the minimal number of 
physical scales to describe these baryons. Recently in a QCD-inspired model of mesons, 
besides the contact interaction, motivated by the hyperfine interaction between the quarks, 
a Coulomb-like potential was considered as well. Going beyond the hyperfine interaction 
itself, the zero range interaction mimics those aspects of QCD that binds the constituent 
quarks in the meson, which in this paper was used to build the baryon. We observed a 
surprising reproduction of the trend and magnitude of the binding energies as a function 
of the distint quark mass. Our conclusion, while unexpected, still cares of more detailed 
analysis which includes the quark spin and the Coulomb-type interaction. The results shown 
here, give a support for the extension of the QCD-inspired model of Refs. |5|U to baryons. 

As the feature we addressed here are valid irrespective of the composite hadron nature, 
we think that our conclusion of the flavor dependence of baryonic masses may still hold in 
a more realistic model. However, we have to stress that the present model is based on the 
notion of constituent quarks, and it only uses the constituent quark mass which in fact, 
because spin effects are averaged, does not really enter into spin dependent interactions in 
our calculations. When spin dependent interactions are present, constituent quark models 



can lead to conflict with data, as has been shown in Ref. |17] for the proton spin measured in 
deep inelastic scattering. Therefore, the extension of our conclusion for the case that spins 
are no longer averaged has to be cautious to avoid conflict with spin data. 
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(Fundagao de Amparo a Pesquisa do Estado de Sao Paulo) and TF thanks FAPESP and 
CNPq (Conselho Nacional de Pesquisa e Desenvolvimento of Brazil). 
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APPENDIX A: TWO-QUARK SCATTERING AMPLITUDE 

The two-body scattering amplitudes r qq (M qq ^j and tq q (Mq^J are the solutions of the 
Bethe-Salpeter equations in the ladder approximation, represented diagrammatically in fig- 
ure 5, for a contact interaction between the quarks [f7|,|I5|. In this case the solution is just 
given by the infinite sum of the product of "bubble" -diagrams (figure 6) multiplied by powers 
of the bare interaction strength. The result is given by the geometrical series: 



1~otq 



(ML) = l —f r , (Al) 



where a = q or Q and A is the bare interaction strength. The function B aq (jM^j is the 
"bubble" -diagram represented in figure 6, 

/ 2 \ r d 4 k i i 

B aq [M aq ) = J (2?r) 4 (fc2 _ m 2 + -) ftp- _ k)2 _ m 2 + ~) > 

where the total four-momentum of the quark pair is P aq with P^ q = M^ q . 

The four-dimensional integration in Eq. (|A2|) is performed in light-front variables. First, 
the virtual propagation of the intermediate quarks is projected at equal light-front times 
f?|,|T2|j, by analytical integration over k~ in the momentum loop. The non-zero contribution 
to Eq. QS^) comes from < k + < P+ q for P+ q > 0, 



i r dk+d 2 k ± [P+ - k+j 6 (fc+) 

Ba g (M aq ) = ——3 / —7—— T7\ — - yr— T- ( A3 ) 



2 (271)" J k+ (P+ - k+) p+ _ kj-rr^ 



Now, we introduce the invariant quantity x = and the relative momentum 

K± = (l- x)kj_ - x(P aq - k)± , 

in Eq. (|A"3|), which gives: 

n (m 2 \ 1 f dxd2R ± 9 ( 1 - x )°( x ) (\ A \ 

aq V- 2 (27T) 3 J X (1 - X) M 2 _ ^i+(m|-mg)x+mg " ^ 

1V1 aq x(l-x) 

The function B aq [Maq) has a log-type divergence in the transverse momentum integra- 
tion. The renormalization T aq (M^ q ) is done by taking into account the physical information 
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of the interacting light-quark pair system, that we suppose has a bound state. Using this 
physical condition to define the two-quark scattering amplitude we have studied the nucleon 



same renormalization condition. 

The pole of the light-quark scattering amplitude, T qq (M qq ) is found when M qq is equal to 
the mass of the bound qq pair, Thus, the bound state pole of the scattering amplitude 
demands that 



which is enough to render finite the scattering amplitudes r qq and tq q . Thus, the bare 
strength of the effective contact interaction between the constituent quarks q and Q does 
not depend on flavor. In this manner, we extend the flavor independence of the gluon 
interaction of the fundamental QCD Lagrangian to the effective interaction. 
The final equation for the two-quark scattering amplitude is: 



in the three-quark light-front model ||. This model fitted, simultaneously, the proton mass, 
the charge radius and the electric form factor below 2 (GeV/c) 2 ||. Here, we just use the 




(A5) 




(A6) 



The log-type divergence of r aq is removed by the subtraction in Eq. (|A6| ) . 
In particular, the analytical form of Eq. ( |A6| ) is |7| , 




(A7) 



for < M qq < 2m q , which is enough for the integration in Eq.(^). 
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TABLES 

TABLE I. Vector (pseudo-scalar) meson and current quark masses from Ref. [16]. The 
estimated values of the current quark masses (last column) are explained in the text. Values 
quoted in GeV. 



Meson 




M M 


P 0) 


ud 


0.768 (0.138) 


K* (K+) 


us 


0.892 (0.494) 


D* (D°) 


Tic 


2.007 (1.865) 


B* (B+) 


lib 


5.325 (5.279) 



quark m" r (Ref. [16]) m c y ,rr 

u,d (1.5-9)xl0- 3 

s 0.06 - 0.17 0.124 

c 1.15 - 1.30 1.115 

b 4 - 4.4 4.557 



TABLE II. Low-lying spin 1/2 baryon experimental masses from Ref. [16] and binding energies 
from Eq. (||). Values quoted in GeV. 



Baryon 



I(J P ) 



M B 



P 
A 

A+ 



uud 
uds 
udc 
udb 



2V2 



0(| + ) 
0(| + ) 

o(D 



0.938 
1.115 
2.285 
5.624 



0.214 
0.161 
0.106 
0.085 
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FIGURES 




FIG. 1. Diagrammatic representation of Eq. (1). The black bubble represents the two-quark 
scattering amplitude. 




FIG. 2. Diagrammatic representation of Eq. (2). The black bubbles represent the two-quark 
scattering amplitudes. 
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FIG. 3. Binding energy (Bjj) as a function of the mass of the respective low-lying hadron 
(Mh)- Experimental data of pseudo-scalar mesons from Table I (full squares). Experimental data 
of the spin 1/2 baryons from Table II (full circles). 
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FIG. 4. Binding energy of the low-lying spin 1/2 baryon states as a function of the respective 
ground state mass of the baryon. The results of the light-front Faddeev model calculation are 
shown by the solid line. Full circles are the data from Table II for the nucleon, A , A+ and A°. 
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+ 



FIG. 5. Bethe-Salpeter equation in ladder approximation for the two quark-scattering ampli- 




T 



tude from Eq. (Al). The contact interaction is represented by the dot. 




FIG. 6. "Bubble" -diagram representing Eq. ([A2|) . 
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